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Abstract. With the advent of Chandra and XMM, X-ray 
lines spectroscopy of AGN outflows is off to an exciting 
start. In this paper we illuminate some of the complica- 
tions involved in extracting the outflow's physical condi- 
tions (ionization equilibrium, total column density, chem- 
ical abundances) from such spectroscopic data. To do so 
we use the example provided by high-quality FUSE and 
HST observations of the outflow seen in NGC 985. 

We show how simple determinations of the column 
density in the UV absorption lines often severely under- 
estimate the real column densities. This is due to strong 
non-black saturation of the absorption troughs, where in 
many cases the UV line profile mainly reflects the velocity- 
dependent covering factor rather than the column density 
distribution. We then show that underestimating individ- 
ual ionic column densities by a factor of 5 can cause a 
two orders of magnitude error in the inferred total column 
density of the outflow. In some case this will be enough 
to associate the UV and X-ray absorber with the same 
outflowing material. 

Finally, we note that the UV spectra of NGC 985 have 
10-20 times the resolution and 2-5 x the S/N per resolu- 
tion element compared with the best available Chandra 
spectra of similar objects. Therefore, in the X-ray band 
non-black saturation and velocity-dependent covering fac- 
tor effects will only become abundantly clear using vastly 
more capable future X-ray telescopes. However, by taking 
into consideration the lessons learned from the UV band, 
we can greatly improve the quality of physical constraints 
we extract from current X-ray data of AGN outflows. 



1. Introduction 

Outflows in Seyfert galaxies are evident by resonance line 
absorption troughs, which are blueshifted with respect to 
the systemic redshift of their emission counterparts. Ve- 
locities of several hundred km s _1 (Crenshaw et al. 1999; 
Kriss et al. 2000) are observed in UV resonance lines 
(e.g., Civ AA1548.20,1550.77, Nv AA1242.80,1238.82, 
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Fig. 1. HST spectrum of NGC 985, taken in 1999 Febru- 
ary 17 with the STIS G140M grating (resolution ~ 
10,000), total exposure 3800 seconds. The absorption 
troughs marked with 1-6 in both Lya and N v arise from 
the outflow emanating from this AGN, Galactic absorp- 
tion is marked with a "G." 

O vi AA1037.62,1031.93 and Lya), as well as in X-ray res- 
onance lines (Kaastraet al. 2000; Kaspi et al. 2000). Simi- 
lar outflows (often with significantly higher velocities) are 
seen in quasars that are the luminous relatives of Seyfert 
galaxies (Weymann et al. 1991; Korista Voit & Weymann 
1993; Arav et al. 2001a). Reliable measurement of the ab- 
sorption column densities in the troughs are crucial for 
determining the ionization equilibrium and abundances of 
the outflows, and the relationship between the UV and the 
ionized X-ray absorbers (Mathur et al. 1995; Elvis 2001). 

In the last few years our group (Arav 1997; Arav et 
al. 1999a; Arav et al. 1999b; de Kool et al. 2001; Arav 
et al. 2001) and others (Barlow 1997, Telfer et al. 1998, 
Churchill et al. 1999, Ganguly et al. 1999) have shown that 
in quasar outflows most lines are saturated even when 
they are not black. We have also shown that in many 
cases the shapes of the troughs are almost entirely due to 
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changes in the line of sight covering as a function of veloc- 
ity, rather than to differences in optical depth (Arav et al. 
1999b; de Kool et al. 2001; Arav et al. 2001a). As a conse- 
quence, the column densities inferred from the depths of 
the troughs are only lower limits. 

These results have led us to suspect that many of the 
reported column densities in Seyfert outflows are likely un- 
derestimated. Recently, we have re-analyzed the HST high 
resolution spectroscopic data of the intrinsic absorber in 
NGC 5548 (Arav, Korista & de Kool 2002) and found that 
the C iv absorption column density in the main trough 
is at least four times larger than previously determined. 
Furthermore, we have shown that similar to the case in 
quasars, the shape of the main trough is almost entirely 
due to changes in covering fraction as a function of veloc- 
ity, and is not due to differences in real optical depth. 

In this paper we analyze the outflow in another 
Seyfert 1 galaxy, NGC 985. Similar to the NGC 5548 case, 
the apparent column density (Lya in this case) is a severe 
underestimate of the real column density. Furthermore, 
the resulting errors in the inferred values of the total hy- 
drogen column density ( N H ) and the ionization parameter 
(U) of the absorbing gas are far larger. For the specific 
case of NGC 985, the real N H and U values may 
allow the X-ray and UV absorbers to arise from 
the same gas. We conclude by highlighting potential dif- 
ficulties expected in the X-ray spectral analysis of AGN 
outflows, based on the UV investigation shown here and 
elsewhere. 

2. The case of NGC 985 

NGC 985 is a bright (V = 13.8 for UV see Fig X-ray: 
< F 2 -iokev >= 2x 10~ n cgs, Nicastro et al. 1998), nearby 
Seyfert 1 galaxy, which shows a strong and complex sys- 
tem of warm absorbers (Njj = 10 22 and 10 23 cm" 2 ) in its 
ROSAT-PSPC (Brandt et al. 1995) and ASCA (Nicastro 
et al. 1998). In figure | we show the available HST data 
of the intrinsic UV absorber in NGC 985, which covers 
the Lya and Nv troughs. To find the physical conditions 
of the absorbing gas we first extract the apparent column 
densities for the Lya and Nv troughs for the outflow, and 
then use these values to compute the N H and U of the 
absorbing gas. 

2.1. Apparent column density extraction and the 
resultant ionization equilibrium determination 

Figure || shows an unsaturated absorption model for the 
Lya trough consisting of an unabsorbed emission model 
convolved with six Gaussian absorption components. The 
fit is very good and the main discrepancy (seen around 
1266 A) is due to our restrictive assumption that each 
kinematic component for Lya and N v will have the same 
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Fig. 2. Lya absorption (histogram) and a simple Gaussian 
absorption model (dotted line) consisting of an emission 
model (solid line) convolved with six Gaussian absorption 
components. To avoid excessive freedom in the model, the 
centroid of each Gaussian lies at the same velocity as the 
corresponding Gaussian we used in modeling the N v ab- 
sorption (see Fig||). 
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Fig. 3. A unsaturated absorption model for the Nv dou- 
blet. We vary the Gaussian parameters to match the ab- 
sorption in the blue doublet component and then scaled 
their depth by half to model the absorption in the red 
doublet component. 

velocity centroid. We note in passing that this assumption 
holds quite well in this case, since four of the Lya com- 
ponents are fitted very well using the velocity centroid 
of the Nv components. For absorption component 4, the 
Gaussian parameters yield log(Njj j) = 14. 

Figure |^ shows an unsaturated absorption model for 
the N v absorption. Components 4 and 5 are not saturated 
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(or only slightly so) since their depths in the red and blue 
doublet lines are consistent with the expected 1:2 optical 
depth ratio (the ratio of their oscillator strengths). For 
absorption component 4, the Gaussian parameters yield 
log(Njv v) = 14.2 

It is usually assumed that the plasma is in photoion- 
ization equilibrium and therefore the relevant physical pa- 
rameters are the total hydrogen column density (N H ) and 
the ionization parameter (U) of the absorbing gas. In or- 
der to determine these we use the photoionization code 
CLOUDY (Fcrland 1996) to try to find combinations of 
these parameters which will reproduce the measured col- 
umn densities. 




-4-2 2 

log(U) 



Fig. 4. Ionization parameter U and total column density 
Njj for component 4 in the NGC 985 outflow, based on the 
apparent column densities extracted from the fits shown 
in figures || and ||. All quantities are shown in log form. 

In figure |] we show photoionization solutions for ab- 
sorption component 4. The presentation method is based 
on the grid models approach developed in Arav et al. 
(2001b), using the same ionizing continuum. We plot 
curves of constant Ni on (determined from the observed 
troughs) in the N H /U plane. That is, within the range of 
plotted parameters, any combination of N H and U that 
falls on the curve yields the desired column density value 
(for an assumed set of elemental abundances, solar in our 
case). The parameters of the absorber are determined by 
the crossing of the different Ni on curves. We note that with 
two given curves we can obtain two solutions, since two 
crossing points are possible. By coincidence, we obtained 
a single solution in this particular case. The inferred N H is 
at least two orders of magnitude too small to account for 



a typical warm absorber column density (Kaastra et al. 
2000; Kaspi et al. 2000). 

2.2. Non-black saturation 
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Fig. 5. Modeling the Ly/3 absorption with the same six 
Gaussians absorption-model used in modeling the Lya 
trough (see Fig. ||). The upper dotted line is based on 
assuming no saturation in Lya. That is, the maximum 
optical depth of each component was multiplied by the 
ratios of oscillator strength and wavelength for Ly/3 and 
Lya, that is by [f ^]i,y/3 / [f ^]i,ya- The poor fit in all compo- 
nents is indicative of saturation in Lya. The lower dotted 
line is an absorption model assuming full Lya saturation 
(i.e., using the same Gaussians with no scaling); moderate 
saturation is evident in components 4, 5 and 6 and small 
saturation is inferred for components 2 and 3. 

However, our ionization equilibrium inferences are cru- 
cially dependent upon the reliability of the Nv and Lya 
column densities inferred from our absorption models 
(Figs. H and ||). The evidence we have collected over the 
last few years suggests that it is impossible in princi- 
ple to extract a reliable column density from one 
singlet line (see § 1 and Arav et al. 2002). It is therefore 
quite possible that our Lya inferred H I column density is 
greatly underestimated. Luckily we possess FUSE spectra 
that covers the Ly/3 and Ly7 outflow troughs in NGC 985, 
and can directly check the reliability of the Lya model. 

Figure [| shows the results of applying the Lya absorp- 
tion model to the Ly/3 trough. Our very good model-fit for 
the Lya trough (Fig. ||) fails completely when we apply 
it to the Ly/3 trough assuming no saturation (upper dot- 
ted line). A model assuming complete saturation (lower 
dotted line) fits the data considerably better, but since 
it over-predicts the depth of the troughs we determine 



4 



Nahum Arav,: AGN Outflows: Analysis of the Absorption Troughs 









i 




\j ' 

1 / 

1 


22 


- 


\ N V 
\14.2 




U = 0.8 J/ / 
N H =21.4/f / 

/ "' 7 


1 

1 

1 


'"a 








/ / 
/ : ' / 


- 


SIT 
O 






/ / 


/ •' / 
/ .'/ 

• ' U=-0.3 
■/ N H =19.3 




18 




H I , 
14.7/ 

i/ i / 


/ / 

/H I 
/ 14.0 

/ 

1 1 


ill 






4 




-2 





6 



log(U) 



Fig. 6. Similar to Fig. |4| showing the ionization equilib- 
rium solutions for both the real (solid line) and apparent 
(dashed line) H I column density. The apparent H I column 
density was extracted from the HST/STIS data alone as- 
suming no saturation, while the real H I column density 
that is fivex larger than the apparent one is based on 
the FUSE data. A factor 5 difference in the H I col- 
umn density increases the inferred total absorber 
column density by more than two order of mag- 
nitudes. We note that the lower crossing-point solution 
(the one without the dot) can be excluded, since we do not 
detect absorption from low ionization species that would 
be present in that case. 

that the Ly/3 is not heavily saturated. This conclusion is 
confirmed by the shape of the Ly7 trough where all the 
components except number 6 are significantly shallower. 
The apparent H I column density of component 4, as mea- 
sured from the Lya trough, is 5x smaller than its real Hi 
column density deduced with the constraints imposed by 
the Ly/3 and Ly7 troughs. In passing, we mention that the 
FUSE spectra also show saturation in the O vi lines. 

3. Realistic ionization equilibrium 

Figure |^ illustrates the consequence of ignoring the possi- 
bility of saturation in Lya and therefore underestimating 
the Hi column density in component 4 by a factor of five. 
The physical solution that uses the real H I column density 
results in an absorber with 100 x higher Nh and roughly 
a 10 x higher ionization parameter than the one using the 
apparent column density. The relative difference between 



the two solutions is a robust key point of this investiga- 
tion. Absolute values of the correct Nh and U are far less 
accurate due to the systematic uncertainties in chemical 
abundances and incident ionizing spectrum, as well as the 
reliance on only to ionic column densities. Nonetheless, we 
point out that he solution based on realistic Hi column 
density is compatible with physical parameters inferred 
for warm absorbers in general (Kaastra et al. 2000; Kaspi 
et al. 2000), while the apparent one is much too low in 
both N ion and U. 

4. Summary and Relevance to analyses of X-ray 
warm absorbers data 

1. Non-black saturation is prevalent in UV troughs of 
AGN outflows, even when the trough is fully resolved 
spectroscopically. This suggests that the same phe- 
nomenon occur in the X-ray troughs of the same ob- 
jects. In addition, a single covering fraction per absorp- 
tion component is a poor approximation, since cover- 
ing fraction is often strongly velocity-dependent. 

2. Errors in column density estimates can be greatly am- 
plified in the inferred ionization equilibrium for the ab- 
sorbing gas. In the case of component 4 in NGC 985, 
correcting the Nh i by a factor of five resulted in a two 
orders of magnitude higher Nh and an order of mag- 
nitude higher U. These differences are large enough 
to make UV component 4 compatible with being the 
same gas as the X-ray warm absorber. 

3. Many line series are covered by the X-ray band, thus 
allowing for excellent saturation diagnostics and thus 
the extraction of reliable column densities. However, 
care must be taken to allow for non-black saturation, 
which is very different from saturation due too spec- 
troscopically unresolved absorption components. 
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